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The Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase, a superconducting state with non-zero total 
momentum Cooper pairs in a strong magnetic field, was predicted more than 50 years ago and now 
becomes an important concept in many branches of physics. However, no unambiguous experimental 
evidences for the existence of FFLO phases have been observed yet. Recently, spin-imbalanced 
ultracold degenerate Fermi gases have emerged as a new powerful platform for the observation of 
FFLO phases due to their high experimental controllability and the lack of disorder. However, 
in three dimensional degenerate Fermi gases where the FFLO phases can be described within the 
simple mean field theory, the parameter region for the FFLO phases is too small to be observed in 
experiments. Recently, we showed that the Rashba type spin-orbit coupling and in-plane Zeeman 
field in three dimensional degenerate Fermi gases provide a more efficient way to create the FFLO 
phase because (1) the parameter region for the FFLO phase is greatly enlarged in the phase diagram 
and (2) the FFLO phase is stabilized due to the enhanced energy difference between FFLO phase 
and conventional Bardeen-Cooper-Schrieffer (BCS) phase. In this work we investigate the FFLO 
phase in two dimensional spin-orbit coupled degenerate Fermi gases using the mean field theory, with 
special concerns on the physical origin of the FFLO state. The basic properties of the FFLO phase 
are discussed and the phase diagram is obtained at zero temperature. The symmetry of different 
quantum phases is examined, which provide important basis for the experimental observation of 
FFLO phases using the time-of-flight imaging. 

PACS numbers: 05.30.Fk, 74.20.Fg, 74.25.Dw 
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I. INTRODUCTION 



In 1964, just shortly after the great success of Bardeen- 
Cooper-Schrieffer (BCS) theory for superconductivity 1 , 
Fulde and Ferrell^, and Larkin and Ovchinnikov 3,4 
(FFLO) independently demonstrated that a new type of 
superconducting state, which is characterized by Cooper 
pairs with nonzero total momentum and spatially mod- 
ulated order parameter, may exist in certain regime of a 
clean superconductor under a strong magnetic field. The 
superconducting state is now known as the FFLO state. 
For the BCS cooper pair—, the pairing takes place be- 
tween electrons with opposite momentum and spin, i.e., 
k j" and — k J,. Therefore when the magnetic field exceeds 
certain critical value, the superconductivity is destroyed 
due to Pauli paramagnetic depairing effect. As a conse- 
quence, the magnetism and superconductivity generally 
cannot coexist for the BCS type superconductor, which 
also fail to explain the origin of magnetic superconduc- 
tors. The story is totally different for the FFLO phase 
because these two different orders naturally coexist, or 
more precisely, the FFLO phase arises from the inter- 
play between magnetism and superconductivity. This 
important feature makes the FFLO phase a central con- 
cept for understanding many exotic phenomena in differ- 
ent physics branches, ranging from unconventional solid 
state superconductors (e.g., layered 5 ^, heavy-fermion£~— , 
organi o 10 ' 11 superconductors, etc.), to chiral quark mat- 
ter in quantum chromodynamics, and to neutron star 
glitches in astrophysic o 12 ' 13 . In the past several decades, 



great efforts have been made to unveil this novel quan- 
tum phase, and a lot of exotic signatures that may be 
related to the FFLO phase have been observed. How- 
ever, until now, the unambiguous experimental evidences 
for the existence of FFLO phases are still lackin g 12 ' 14 . 
There are several reasons for that: firstly, the existence 
of FFLO phase requires very stringent conditions; and 
secondly, the direct probing of periodic oscillation of the 
order parameter is challenging; and thirdly, the disorder 
effects in the superconductor induce strong scattering be- 
tween different momenta that destroys the superconduct- 
ing pairing 
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The experimental advances of population imbalanced 
ultracold Fermi gases have the potential to elucidate this 
long sought problem. The ultracold atomic system pos- 
sesses some remarkable advantages over the counterpart 
in solid state systems due to its high controllability and 
tunabiUtyil-iS. This new platform in principle can be 
made disorder free, and the disorder, if necessary, can be 
introduced to the system in a controllable manner— ~—. 
On the experimental side, the superfluidity of the Fermi 
gas can be characterized by the generation of vortices 
when the gas is rotated 2 ^., and the momentum of the 
Cooper pair in the FFLO phase can be directly probed 
using the time-of-flight imaging 24,25 , while in solid state 
the direct observation of the FFLO phase is challeng- 
ing. The main problem is that the FFLO phase only 
exists in an extremely narrow parameter regime in three 
dimensional degenerate Fermi gase o 24 ' 26 , therefore in ex- 
periments the FFLO phase is generally missed out. For 
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FIG. 1: (Color online) Basic physical picture for the emer- 
gence of FFLO phase in ultracold atomic system, (a) shows 
the creation of FFLO phase with Zeeman field, which is equiv- 
alent to the physics by controlling the population imbalance. 
The mismatched Fermi surface makes the pairing with oppo- 
site momentum and spin challenging, thus in some parameter 
regime the FFLO phase becomes energetically favorable. This 
mechanism is inefficient because the Pauli paramagnetic de- 
pairing effect under strong Zeeman field destroys the Cooper 
pairs, thus the FFLO phase only survives in a very narrow 
parameter regime, (b) shows a new mechanism for the gen- 
eration of FFLO phase using the deformation of the Fermi 
surface. Such deformation of the Fermi surface can be con- 
structed by an in-plane Zeeman field and SO coupling. The 
center of the Fermi surface is also shifted due to the Zee- 
man field, and as a consequence, the BCS type of pairing al- 
ways becomes challenging even with a small Zeeman field. A 
very large FFLO phase can be found in the parameter space. 
The deformation of the Fermi surface breaks the rotational 
symmetry of the system, thus has the potential to create the 
FFLO phase with a unique FFLO momentum Q. 



instance, in recent experiments with population imbal- 
anced Fermi gases by Zwierlein et a / 27 ' 28 only the phase 
transition from BCS superfluid to normal gas has been 
observed, while in another parallel experiment by Par- 
tridge et a& the phase separation phase, which is also 
known as the breached pair—, has been observed. There 
are two basic reasons why the FFLO phase is missed out 
in these experiments. Firstly, the parameter space for 
the FFLO phase is extremely narrow; and secondly, the 
energy difference between the FFLO ground state and 



(a) Without SOC (b)WithSOC 




FIG. 2: (Color online) Influence of SO coupling on the emer- 
gence of FFLO phase. In (a) , the FFLO phase emerges when 
the Zeeman field exceeds certain critical value. However, the 
FFLO ground state energy only depends on the magnitude of 
the FFLO momentum Q but is independent of the direction, 
thus in principle all the momenta with the same magnitude 
(|Q| = const, see the circle in (a)) are the system's ground 
state. In the SO coupled system Q becomes unique in (b) 
because the SO coupling and in-plane Zeeman field break the 
rotational symmetry of the model. 



BCS excited state is much smaller than the tempera- 
ture. Therefore even the parameters for the FFLO state 
have been reached, the FFLO ground state cannot be 
observed. 

In this paper, we show that the above two limitations 
can be overcome using spin-orbit (SO) coupled degen- 
erate Fermi gases with an in-plane Zeeman field. Dif- 
ferent from SO coupling in solid state systems which 
arises from spatial inversion symmetry breaking^, SO 
coupling in ultracold atoms is generated through laser 
couplings^ 2 - - — between hyperfine states, therefore differ- 
ent types of SO coupling can in principle be created by 
carefully choosing different laser configurations. Experi- 
mentally, certain type of SO coupling has been realized 
using Raman coupling between hyperfine states for both 
Bose and Fermi gases^ 8 . - — . In this paper we use Rashba 
type of SO coupling as an example to demonstrate our 
basic ideas. The Zeeman field is assumed to be along 
the x direction, which is in the SO coupling plane. The 
basic idea for the emergence of the FFLO phase is illus- 
trated in Fig. QJ>, which also applies to other types of 
SO coupling. For the Fermi gas with only Zeeman field, 
see Fig. the two mismatched Fermi surfaces always 
form concentric circles, therefore for the s-wave pairing, 
the up- and down-spins acquire different Fermi momen- 
tum, i.e., k + Q/2,ft and -k + Q/2 J|, with ft and JJ- 
spins in the pseudospin representation and Q as the to- 
tal momentum of the Cooper pairs. The free energy of 
the system satisfies the rotational symmetry 

F(Q) = F(|Q|). (1) 

This is because the Zeeman field only fixes the direction 
of the spin, but does not fix the direction of the mo- 
mentum axis, therefore the free energy should have the 
rotation symmetry for the momentum. It can also be un- 
derstood from the fact that the total free energy depends 
on k 2 , Q 2 and k ■ Q, thus the summation over k should 
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be independent of the direction of Q, see Ref. 24 Phys- 
ically, Eq. ([T]) means that the total momentum of the 
Cooper pair can take any direction therefore the ground 
state FFLO phase is infinity-fold degenerate. Generally 
in the numerical simulation, we artificially set Q along 
a particular direction and demonstrate that the FFLO 
phase indeed has a lower energy than the regular BCS 
superfluid (Q = 0). Due to the Pauli paramagnetic de- 
pairing effect, the FFLO phase only survives in a very 
narrow parameter regime, see also the numerical results 
in Fig. ilk. 

The physical picture is totally different when a Rashba 
SO coupling is included, as schematically shown in Fig. 
[Tb. In this case the Fermi surface is deformed and the 
center of the Fermi surface is no longer located at k = 0. 
Notice that the deformation of the Fermi surface depends 
strongly on the direction of the SO coupling and Zeeman 
field. For the model we consider here, the deformation 
is along the y direction. In the pseudospin representa- 
tion (the eigenstates of single particle Hamiltonian) , we 
have both singlet pairing and triplet pairing. The triplet 
pairing is almost intact by the Zeeman field, thus the 
FFLO phase can be observed in a much larger parameter 
regime. The deformation of the Fermi surface makes the 
FFLO phase always energetically favorable even with a 
small Zeeman field. We have numerically confirmed that 
the FFLO momentum Q should be along the y direction, 
and the minimization of the total free energy should give 
a single Q value, therefore the FFLO ground state should 
be unique, see also Fig. \Bp for the movement of Q in the 
Q — A space. For the symmetry argument, see Appendix 
A. 

Generally, the mismatch of the Fermi surface is the ba- 
sic route to the FFLO phase, and such mismatched Fermi 
surface can be created by population imbalance^—, Zee- 
man fields or mass imbalanc o 44 ' 45 . In this work, to- 
gether with our previous work 4 ^, we demonstrate that 
the FFLO phase can be created more efficiently through 
the deformation of the Fermi surface, which can be 
constructed by SO coupling, or, non-Abelian gauge 
field^-^L 4 !, and Zeeman field. Notice that the gener- 
ation of non-Abelian gauge fields is a subject of intensive 
investigations in ultracold atoms in the past decade, see 
a recent review in Ref. |U For this new route, the Zee- 
man field is still needed, otherwise, the system has time- 
reversal symmetry and the band structure should satisfy 
E>kf = i?_k,i, which means that two atoms with oppo- 
site momentum and spin on the Fermi surface can always 
form BCS Cooper pairs efficiently, leading to BCS super- 
fluid, instead of FFLO phases. Our route here, however, 
shows that the FFLO phase may become energetically 
favorable with a small Zeeman field. 

In Ref. Hi| we have considered a three dimensional de- 
generate Fermi gas with Rashba SO coupling and an in- 
plane (in the SO plane) Zeeman field. The FFLO phase 
cannot be observed in three dimensional solid state sys- 
tems because the orbital effect is unavoidable no matter 
which direction the external magnetic field is applied. 



The low dimensional systems in solid state materials pro- 
vide the opportunity for observing FFLO states because 
the orbital effect can be suppressed by applying a mag- 
netic field parallel to the system, which only leads to a 
Zeeman field. In ultracold atomic systems, the Zeeman 
field is generated through different methods without in- 
volving the orbital effect, therefore the FFLO state may 
be observed even in three dimension. 

Here we want to emphasize the difference between our 
proposal and the results in solid state systems 4 ^. In 
solid state materials, it has been found that the Rashba 
type SO coupling in superconducting materials without 
inversion symmetry can enhance the FFLO phase 4 ^ by 
increasing H C 2, the so-called upper critical field of the 
type-II superconductor—. In solid state materials, the 
SO coupling energy aKp and Zeeman splitting generally 
are much smaller than Ep, therefore the deformation of 
the Fermi surface is very small, and cannot be the main 
driving force for the FFLO phase. In our model the SO 
coupling energy, Zeeman field and Ep are generally of 
the same order, thus the deformation of Fermi surface, 
as shown schematically in Fig. [TJd, is significant and pro- 
vides a major driving force for the formation of the FFLO 
phase. Moreover, in solid state materials, the disorder ef- 
fect is one of the major obstacles for the observation of 
FFLO phased, while the ultracold atomic systems can 
be made disorder free. 

The rest of this work is organized as following. We 
present our mean field treatment of the SO coupled de- 
generate Fermi gas with an in-plane Zeeman field in sec. 
HT1 Then we plot the phase diagram and discuss its basic 
properties in sec. Mil We discuss the problems related 
to the measurement of the FFLO phase in sec. IIVI At 
last we conclude in sec. |Vj The symmetry of the effective 
Hamiltonian is discussed in Appendix A, and the treat- 
ment of the same model using the quantum field theory 
is presented in Appendix B. 

II. PHYSICAL MODEL, HAMILTONIAN, AND 
GAP EQUATION 

We consider a two dimensional degenerate Fermi gas 
with Rashba-type SO coupling and an in-plane Zeeman 
field. The two dimensional degenerate Fermi gases can 
be constructed by applying a strong standing wave along 
the third direction, and have been realized in recent 
experiments^. The two dimensional SO coupled Fermi 
gases can be described as 

H = ^ 4..g [£ k <7 +a(k x a y -ky<j x )- ha x ] c k , ff ' +Vmt , (2) 

kcrcr' 

where a is the SO coupling strength, a x and a v are the 

Pauli operators, £ko- = ^ — Mi an d k = (k Xl k y ). The last 
term corresponds to the s-wave scattering interaction 

V int =9 c Lt C P2 4 C P3 4 c P4,-b (3) 

Pl+P2=P3+P4 
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where pi + p 2 = P3 + P4 due to the conservation of 
the total momentum during the scatting process, g is the 
scattering interaction strength. 

When atoms form Cooper pairs with a finite total 
momentum Q, the scattering process in Eq. [3] can 
be simplified with pi = k + Q/2, P2 = k + Q/2, 
P3 = P + Q/2, and p 4 = p + Q/2. When Q = 0, the 
Cooper pairs are formed between two atoms with op- 
posite momentum and spin, recovering the conventional 
BCS type of pairing. Denote /3 p = gCp + Q/ 2 ,ic_ p+ Q/ 2: f, 
the interaction term can be rewritten as V; ni = 

t t 

9 2^k,p C k+Q/2,f C -k+Q/2,| C P+Q/24 C -p+Q/2,t 



Ek,p 



The interaction term can be decoupled 



using the standard mean-field method 




(4) 



FIG. 3: Particle-hole symmetry in our model. The particle- 
hole symmetry establishes the direct relation between ?/>k for 
Hefi(k) with eigenvalues _Bk to E^k for £f e ff(— k) with eigen- 
values —-Ek- Due to the lack of chiral symmetry^, the eigen- 
values of f/cff (k) may not appear with pairs. The particle-hole 
correspondence is shown via dashed line. 



where the order parameter in the momentum space A = 
Ep9( c p+Q/2,l c -P+Q/2,t>- The interaction term now re- 
duces to 



^A/?£ + A*/3 k 



IAI 



(5) 



Notice that c pt7 = Jdxc cr (x)e lpx , then we have the 
paring in the real space, 



(c t (x)c;(y)) = J dpdp' (.gc p , t c pU ) 
= A5(x-y)e jQ ' x . 



i(p-x+p'-y) 



(0) 



We see the order parameter defined in the momentum 
space corresponds to a non-uniform pairing in the real 
space. The 5-function arises from the contact interaction. 
The above pairing breaks the time-reversal symmetry, 
which is in consistent with our model because a Zeeman 
field is applied. 

Using the standard Bogliubov transformation, the 
Hamiltonian can be rewritten as 



where the effective Hamiltonian reads as, 



rr _(KQi)) 

v At/ 2x2 



A/ 2x2 
-a y K*{-k)a v 



with 



R(k) = 
diag(l, I) 



/C(k) 



£k+Q/2,f 

h - aR(k)* 



h - aR(k) 
Ck+Q/2,4. 



(7) 



(8) 



(9) 



(k + Q/2) x + i(k + Q/2)„ and 7 2x2 
The basis defined in Eq. [7] is ipKQ 

J 



(c k +Q/2,T> c k+Q/24' C -k+Q/2,i'~ C -k+Q/2,t) ■ The mi - 

nus sign in the last term in the basis is essential to achieve 
the AJ 2x2 type of off-diagonal term in Eq. [5] 



The symmetry of the system is essential for under- 
standing the properties of different quantum phases, and 
directly related to the measurement of the FFLO state 
using the time-of- flight measurement. We first consider 
the particle-hole symmetry of the model. The particle- 
hole operator can be written as £ = AK, where A = <y y T y 
(a y is the Pauli spin matrix and t v is the Nambu particle- 
hole matrix), and K represents the complex conjugate 
operator. A is a unitary operator, and A = A -1 . K is an 
antiunitary operator because K 2 = 1, which is used to 
ensure that £ is an antiunitary operator. Note that the 
system does not have the time-reversal symmetry but has 
the particle-hole symmetry with E 2 = AKAK = AA* = 
A 2 = 1, therefore the system belongs to the symmetry 
class D, see Ref. USE We have 



EHefftk)!]- 1 = AH * ff (k)A 
= A f£*(k)) 



At/ 2x2 
\ v A/ 2x2 -a y K(-k)a y 

—i(Ty\ ( tA^CTy 

ia v 







— ia % 

O y 



-iJC* (k))a y 
— zAcr,, 



K -ia v 1C(-k)) 
-/C(-k)) -A/ 2x2 ' 

-A t / 2X2 (TylC* (k)(Ty / 

= -H cS (-k). (10) 

Here £ establishes a one-to-one correspondence between 
k and — k, therefore if V'k = ( u (k), f(k)) T is an eigen- 
vector of the Hamiltonian H c g(k) with energy E^, then 
ip' k = S-^k = (a y v*(k), — a y u*(k)) T is the eigenvec- 
tor of the Hamiltonian H e s(— k) with energy —E(k), as 
schematically shown in Fig. [3] The particle-hole symme- 
try does not automatically ensure that the eigenenergies 
appear in pairs (E, —E) because the system lacks the 
chiral symmetry^. This is a direct consequence of the 
rotational symmetry breaking for our model due to the 
in-plane Zeeman field. However, for an out-of-plane Zee- 
man field, we observe that all eigenvalues still appear 
with pairs, see more details in Appendix A. For the BCS 
superfluid, we find tr(if e fj(k)) = ^ A E\ — 0, while for 
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the FFLO phase, we find tr( J tf eff (k)) = ^ ^ o. There- 
fore for the FFLO phase, the eigenvalues never appear 
with pairs no matter which direction the Zeeman field 
takes. 

We are interested in the particular case when Q is 
along the y-axis (i.e., Q x = 0). In this case, the sys- 
tem has a discrete Z2 symmetry, i.e., the band struc- 
ture is invariant under transformation k x — > — k x , hence 
E\(k x ,k y ) — E\(—k x ,k y ). Such symmetry can be di- 
rectly probed using the timc-of-flight imaging, see below. 
We also find that a topological phase transition may take 
place when h = ^/ + A 2 for the BCS superfluid, which 
can be obtained following the same procedure in Ref. 
see also Appendix A. In our numerical simulation, we 
find that the Zeeman field and SO coupling generally re- 
duce the order parameter strength, and the superfluid 
is always destroyed before entering the topological BCS 
superfluid regime. The topological superfluid may be ob- 
served when the out-of-plane and in-plane Zeeman fields 
are presented simultaneously. However, since the particle 
and hole may not appear with pairs due to the rotational 
symmetry breaking, it may not be suitable for the ob- 
servation the associated topological excitations, such as 
the zero-energy Majorana Fermion. With only out-of- 
plane Zeeman field, the topological phase transition has 
been intensively discussed in the literature recently^— . 
The model with in-plane Zeeman field generally favors 
the FFLO phase and suppresses the topological super- 
fluid, therefore only regular BCS superfluid phase, FFLO 
phase, and normal gas phase can be observed in our 
model. 

The thermodynamical potential at zero temperature 
reads as 

« = ~ + \ + \ E M-Ex), (11) 
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k.A 



where the Heaviside step function 



Q(x) 




(12) 



E\, A =1, 2, 3 and 4, are the eigenvalues of the effec- 
tive Hamiltonian iJ e ff, whose exact expressions are too 
complex to be presented here. To facilitate the further 
investigation of the FFLO phase, especially the physics at 
finite temperature where the fluctuation effect becomes 
important, we also derive the same thermodynamical 
potential using the quantum field theory, see Appendix 
B. We therefore see the standard mean- field decoupling 
used in Eq. 2] actually corresponds to the Hubbard- 
Stratanovich transformation in the quantum field theory. 
We use the mean field theory as the main theoretical tool 
in this work because it provides a more transparent in- 
terpretation of the FFLO physics. 

The effective scattering interaction g in Eq. [11] in a 
two dimensional Fermi gas should be regularized using 60 



1 



E 



1 



k 2 /m + £ 6 ' 



(13) 
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FIG. 4: (Color online). Phase diagram of the two dimensional 
degenerate Fermi gas in the presence of Rashba-type SO cou- 
pling and an in-plane Zeeman field, (a) results with vanishing 
SO coupling, (b), (c) correspond to the results with SO cou- 
pling ciKf = 0.5 and <xKf = 1.0, respectively, (d) shows 
the phase diagram in the h — oKf plane at Eb = 0.4_Ef- hi 
(h 2 ) defines the boundary between BCS (FFLO) and FFLO 
(normal gas) phase. In our calculation, we do not to observe 
topological superfluid phases. 



where the binding energy Et can be tuned by vary- 
ing the s-wave scattering length through Feshbach 

1 7-1 9 

resonance^—. 

The order parameter A, chemical potential /i, and the 
FFLO momentum Q should be solved self-consistently 
due to the conservation of atom number, i.e., 



= ~n, 



on 

dA 



= 0, 



on 

dQ 



= 0. 



(14) 



Here A and Q are used to minimize the thermodynamical 
potential CI. We consider three different quantum phases: 
the normal phase with A = and Q = (for the normal 
gas Q does not enter the effective free energy, thus can be 
any value. We force Q = 0); The BCS type of superfluid 
with Q = but A ^ 0; and the FFLO phase with Q ^ 
and A ^ 0. When we fix Q = then only BCS type 
of superfluid phase and the normal gas can be obtained. 
Throughout this work, we use two different strategies to 
check the influence of Q on the formation of the FFLO 
phase. In the first strategy, we enforce Q = while in 
the other strategy, we let Q as a free parameter. For the 
results at Q ^ 0, these two strategies yield the energy 
difference between the FFLO ground state and the pos- 
sible BCS superfluid excited state, which is crucial for 
the stability of the FFLO phase at finite temperature. 

Because the Zeeman field is applied along the x- 
axis, the population imbalance should be defined us- 
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FIG. 6: (Color online). Population imbalance P = (o x )/n as 
a function of SO coupling strength. E b = OAEf- 



FIG. 5: (Color online). Evolution of chemical potential (a), 
order parameter (b), and FFLO vector Q (c) as a function of 
the binding energy. In (d) we plot (Jpflo — FBCs)/nEF vs. 
Ei,. In all calculations we set h = 0.8Ef- 

ing the eigenstates of a x , instead of <r z . Since (a x ) — 
Sk( c iL t Ck ^ ~^ n - c )' we therefore have 

P = — =^^^A^ 2X2 -a x ^ 

k,A v 7 

= ^E<Q,a(qI: 2 o^j) ^,Q,a.(15) 

k v 7 

Here V'k.Q.A is the eigenstate of the effective Hamiltonian 
H cS , i.e., -ff c fr0k,Q,A = E\ip^Q^. 

In our calculation, we choose the energy unit as the 
Fermi energy E F of the system without interaction, Zee- 
man field and SO coupling. The corresponding length 
scale Kp 1 is defined through the Fermi momentum K F . 
At finite temperature, the two dimensional system does 
not have the long-range order due to the phase fluctua- 
tion and the relevant physics is the Kosterlitz-Thouless 
transition^. In this paper, we restrict to the physics 
at zero temperature, where the mean-field theory is still 
valid. We have numerically verify that the FFLO phase 
has a smaller free energy when Q//y, as compared to 
that along the x direction, therefore we consider only 
Q = (0, Q). Similar conclusion can be found in three di- 
mensional SO coupled degenerate Fermi gas, see Ref. |H, 
and in two dimensional superconductors where Q 1 B for 
FFLO phases. We notice that the direction of the FFLO 
vector Q is also consistent with the results in solid state 
systems with weak SO coupling, see Ref. HH. 

III. PHASE DIAGRAM 

We first present the phase diagram with different SO 
coupling strength and Zeeman field in Fig. SJ Without 



SO coupling, see Fig. |4jt, we see that the FFLO phase 
only exists in an extremely narrow parameter regime. 
When E b > Q.7Ep, the FFLO phase disappears, thus 
such a phase only can be observed in the weak binding 
energy regime, for instance, E b e (0.15, 0.7)E F . Simi- 
larly, the FFLO phase can also be observed in the three 
dimensional system, see Ref. |46|. However, the FFLO 
phase in three dimensional Fermi gases can only be ob- 
served near the unitary regime within a small parameter 
region, and the small FFLO regime can be easily missed 
out in realistic experiments, which is also one of the main 
reasons why the FFLO phases cannot be observed in re- 
cent experiments in three dimensional Fermi gases2£~— . 
With an increasing SO coupling strength, see Fig. |4Jd for 
aK F = 0.5E F and Fig. ifc for aK F = 1.0E F , we find 
that the FFLO phase regime is greatly enlarged. In the 
strong SO regime in Fig. |4j;, we even observe that the 
phase diagram is almost fully filled by the FFLO phase, 
while the BCS superfluid phase is greatly suppressed and 
only survives in a very small regime. To see the impact of 
SO coupling more clearly, we plot in Fig. |3Jl the phase di- 
agram in the h — aK F plane with E b = 0AE F . We define 
the boundary between BCS superfluid and FFLO phase 
as hi and the boundary between FFLO phase and nor- 
mal gas as /12 for convenience, see Fig. We observe 
hi decreases while increases with the increasing SO 
coupling strength, therefore the FFLO phase is greatly 
enlarged in the strong SO coupling regime. It should be 
noticed that in three dimensional Fermi gases slightly 
decreases with the increasing SO coupling strength 46 . In 
the strong SO coupling region, hi becomes very small, 
but never becomes zero because the the Zeeman field is 
essential for the FFLO phase, which breaks the time- 
reversal symmetry. 

We plot the evolution of chemical potential, order pa- 
rameter and Q as a function of binding energy in Fig. 
[3 where the Zeeman field is fixed to h = 0.8E F . As 
we decreases the binding energy, we observe a sudden 
drop of the order parameter in Fig. [SJd at zero SO cou- 
pling strength due to the Pauli paramagnetic depairing 
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FIG. 7: (Color online). Influence of SO coupling and in-plane Zeeman field on the free energy per particle, F/(uEf), in the 
Q — A plane. The cross symbol in each panel corresponds to the self-consistent solution of Eq. I14I In (a3) and (b3) only 
normal gas can be observed, thus A = 0, and Q can be any value because it doesn't enter the free energy. 



effect, following which there is a small regime that sup- 
ports FFLO phase, see also the solid line in Fig. [5};, 
Q^0. With the increasing SO coupling strength, we see 
that the change of A becomes a smooth function of Eb, 
and in a much larger parameter regime we can observe 
the FFLO phase with non-zero Q. The decrease of the 
order parameter in the FFLO phase regime tends to re- 
duce the free energy, therefore in Fig. [2]the arrows always 
move toward the direction with smaller order parameters 
with the generation of the FFLO phase. 
We further consider the free energy 

F = fl + nn, (16) 

where f2 is the thermodynamical potential defined in 
Eq. [Til regularized by Eq. [13j We denote Ffflo 
as the free energy obtained by letting Q as a free pa- 
rameter, while -Fbcs as the free energy by enforcing 
Q = 0. In the FFLO phase regime, Fbcs represents 
the free energy of BCS excited states, therefore the en- 
ergy difference per particle between Ffflo and Fbcs, 
i.e., 5F = (Ffflo — Fbcs)/iT'Ef characterizes the sta- 
bility of the FFLO phase (i.e., the larger \SF\, the more 



stable FFLO phase). The numerical results are presented 
in Fig. [5]i, where we clearly see the enhancement of \5F\ 
due to the SO coupling. Note that in two dimensional 
Fermi gases the enhanced factor is about 3, which is 
much smaller than that in three dimensional Fermi gases 
(~ 100). 

In Fig. [5p, we see that in the BCS superfiuid regime 
(Eb > 0.7Ep), the order parameter decreases with the in- 
creasing SO coupling, which is in sharp contrast to that 
for SO coupled BEC-BCS crossover with Z direction Zee- 
man field. Generally, with the Z direction Zeeman field, 
the SO coupling plays the role of increasing the density 
of states near the Fermi surface, which increases the or- 
der parameters as well as the critical temperature. With 
an in-plane Zeeman field, the SO coupling plays a to- 
tally different role. To see this more clearly, we plot the 
population imbalance P = (cr x )/n as a function of SO 
coupling in Fig. [51 We see that the SO coupling gener- 
ally increases the population imbalance, thus renders the 
decrease of the order parameter as observed in Fig. [SJd, 
which also hinders the appearance of topological super- 
fluid. In the FFLO phase regime, the order parameter 
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FIG. 8: (color online). Eigenvalues E x (A = 1, 2, 3, 4) of the 
SO coupled degenerate Fermi gas. (a), (b) correspond to the 
typical eigenvalues E\ for the BCS superfiuid with parame- 
ters h = 0.2, Eb = 0.4, akF = 1.0. (c), (d) correspond to the 
typical eigenvalues E\ for the FFLO superfiuid with param- 
eters h = 0.4, Eb — 0.4, akF = 1.0. (e), (f) correspond to 
the typical eigenvalues E\ for the normal gas with parameters 
h — 1.0, Et — 0.4, akF = 1.0. The first column shows the 
results along the x direction, while the second column shows 
the results along the y direction. 



FIG. 10: (color online). Momentum distributions n CT (fe) = 
(cl(k)ccr(k)) and n = n_i + rif for different quantum phases. 
Other parameters are exactly the same as that in Fig. [8] 
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FIG. 9: Gapless excitations for the FFLO phase. Near 
Z?k,A = 0, the dispersion show clear linear dispersion. In two 
dimensional system the linear dispersion is essential to make 
the FFLO phase robust against low-energy fluctuations. 



increases with the increasing SO coupling strength due 
to the formation of the FFLO phase. 

The basic physical picture in Fig. [Tp represents a to- 
tally different and more efficient route to find the FFLO 
state. However, the physical picture along is insufficient 
to explain the emergence of the FFLO phase. To under- 
stand the results more clearly, we plot the free energy per 



particle in the Q — A plane, where the global minimum 
of the free energy marked by the cross symbol in each 
panel corresponds to the self-consistent solution of Eq. 
[141 The left to right shows the influence of the SO cou- 
pling on the formation of the FFLO phase, while the top 
to down shows the influence of the Zeeman field. Note 
that without SO coupling, the free energy is a symmet- 
ric function of Q, therefore in Fig. [7^,2 there are two 
degenerate FFLO ground states at ±Q because we have 
assumed Q is along the y direction in our numerical sim- 
ulation. This is not true in the real system because the 
free energy only depends on the magnitude of \Q\, but 
is independent of the direction of Q due to the rotation 
symmetry, see Eq. [1] The free energy is still a symmetric 
function with respect to Q when h = 0. However it be- 
comes an asymmetric function when both Zeeman field 
and SO coupling strength become non-zero, thus only 
one global minimum can be found in the Q — A plane 
in Figs. [7p2, c2, c3, and the ground state is unique. 
The plot of the free energy in the Q — A plane reflects 
the effect of the symmetry of the effective Hamiltonian, 
see also Appendix A. With the increasing Zeeman field, 
the order parameter decreases due to the formation of 
the FFLO phase, which is in consistent with the basic 
physical picture in Fig. [2] Due to the increase of I12, 
the boundary between FFLO phase and normal gas, we 
observe the FFLO phase in the strong Zeeman field and 
strong SO coupling region in Fig. 0:3. 
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IV. MEASUREMENT OF THE FFLO PHASE 

The three different phases have different properties 
which can be used for the identification of these phases. 
In Fig. [HI we plot the typical band structures E\, A = 1, 
2, 3, 4, for the BCS superfluid, the FFLO phase and the 
normal gas. Due to the rotational symmetry breaking, 
we have to plot the dispersions along the k x and k y axes, 
respectively. For a typical BCS superfluid (Q = 0) in Fig. 
EH and Fig. [5Jd, we see that the system is always gapped 
and the band structure is always symmetric about k = 
for the dispersion along k. x . While along the k y axis, 
such symmetry is absent. For the BCS superfluid we can 
prove that when h 2 < A 2 and A ^ 0, the BCS super- 
fluid is always gapped, see Appendix A. However for the 
FFLO phase, the superfluid becomes gapless along both 
k x and k y axes. Along the k x axis the band structure is 
symmetric about k = 0, but along the k y direction such 
symmetry is broken. For the FFLO phase we observe 
J2\ E\ ^ because Q ^ 0, which is consistent with our 
symmetry analysis in Appendix A. In fact, gapless exci- 
tations are a typical feature of the FFLO phase, as have 
been pointed out in literatur e) 24 ' 61 . In the vicinity of the 
gapless excitation, see Fig. [9j the dispersion becomes 
linear which is essential to ensure that the FFLO phase 
is robust against the low-energy fluctuations. Here we 
should emphasize that not all FFLO phases are gapless. 
The FFLO state may become gapless only when Q is rel- 
atively large, while for a small Q (near the boundary be- 
tween FFLO and BCS superfluid) the FFLO phase is still 
gapped, similar to that in the BCS superfluid. For the 
normal gas the band structure also shows strong asym- 
metry along the k y axis, as seen in Fig. [8^ and Fig. [8f. 

The corresponding momentum distributions n a — 
(c ka Cka) and n — n-f + nj. provide an important tool to 
detect the properties of the FFLO state because they can 
be directly measured via free expansion of the cloud. We 
plot the momentum distributions in Fig. [10] for the three 
different phases presented in Fig. [5] at zero temperature. 
The dispersion properties of the band structures can be 
directly reflected on the corresponding momentum distri- 
butions. We see that for three different quantum phases, 
the momentum distributions are always symmetric about 
k = along the k x direction, while show strong asym- 
metric along the k y direction. However, the sum of the 
momentum distributions n for spin up and spin down 
components still shows perfect symmetry about k = 
along both k x and k y directions. Therefore detecting the 
asymmetry of the superfluid is not sufficient for the iden- 
tification of the FFLO phase. To identify the superfluid 
nature of the FFLO phase, we have to rotate the sample 
to create vortices, which is a direct evidence of superflu- 
idity. Near the boundary between difference phases, the 
fluctuation effect may become significant thus the phase 
boundary region is not suitable for the observation of vor- 
tices. With the large FFLO phase region in our model we 
can safely choose some parameters in the middle of the 
FFLO phase where the fluctuation effect should be min- 



imized. The large FFLO superfluid phase ensures that it 
will not be missed out in future realistic experiments. 

The properties of the FFLO phase may be mea- 
sured using a number of methods developed in ultra- 
cold atom systems, for instance, shot-noise correlation 62 
and density-density correlation measuremen t 63 ' 64 , which 
show a peak at the Cooper pair momentum Q. After 
released from a trapping potential, the free expansion 
of the Fermi cloud has a peak at r = HQt/m, there- 
fore the direct measurement of the FFLO momentum 
Q is possible^. In our model when Q is unique, re- 
peated measurement to determine the FFLO momentum 
becomes possible. In the FFLO phase without SO cou- 
pling, the ground states is independent of the direction 
of Q, then such measurement may become difficult. The 
FFLO phase can also be measured using the Fourier sam- 
pling of time-of- flight images proposed by Duar* 6 - 6 . The 
gapless excitations in the FFLO phase may be observed 
using the Bragg spectroscopy 6 ^. 



V. CONCLUSION 

To summarize, in this paper, we study the possible 
FFLO phase in SO coupled degenerate Fermi gases with 
in-plane Zeeman fields. We show that the FFLO phase 
can be greatly enlarged due to the deformation of the 
Fermi surface. The emergence of the FFLO phase is un- 
derstood from different angles. The properties of the 
BCS superfluid, FFLO phase and normal gas have also 
been discussed and their measurement through the time- 
of-flight imaging is presented. The deformation of the 
Fermi surface provides a more efficient method to gener- 
ate the FFLO phase. On the experimental side, the SO 
coupling has been realized in Bose^r— and Ferm i 42 i 43 
gases, while the later is more relevant to the topics in 
this work. Moreover, the in-plane Zeeman field can be 
naturally created 3 ^ 40 ' 42 ' 43 and controlled at will, thus we 
expect the idea in this work may provide a path for elu- 
cidating the long-standing problem about FFLO phases 
in experiments in the near future. 
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Appendix A: Symmetry of the effective Hamiltonian 

The symmetry of the effective Hamiltonian is one of the 
major concerns in this work. In the following, we discuss 
the symmetry difference between the out-of-plane and 
in-plane Zeeman field models. We consider the following 
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generalized effective model, 



= / AC(k)) A/ 2x2 

yA t / 2x 2 — C7yK*(—V)(T v 



with 



JC(k) 



Ck+Q/2,t 
/l - ai?(k)* 



afl(k)\ 

fQ/24 / 



(Al) 



(A2) 



where £ko- = ^7 — fi + crT. We have applied Zeeman fields 
along Z direction and X direction simultaneously. 

Here we only analyze the symmetry of the BCS super- 
fluid in the generalized model in Eq. IA11 which is more 
relevant to topological excitations such as zero-energy 
Majorana Fermions in vortex cores of the superfluid. We 
hence set Q = 0. For the case with only out-of-plane 
Zeeman field, we set h = 0. The SO coupling term can 
be written as V so = a(k x a y — k y a x ), therefore we define 



The topological properties of the BCS superfluid can 
be calculated from the PfafBan of the following matrix 



W(k) = H cS (k)A = H cS (k) 







-IG, 



(A7) 



Notice that the particle-hole symmetry A_ff* ff (k)A = 

-Feff(-k), where A = A T = At, thus H* s (k)A = 

-AJJrfK-k), and At#J ff (k) = -ff e f ff (-k)At = 
-# eff (-k)A+. We have 

W T (k) = A T iJ c T ff (k) = A+ffeffC-kJAt = -W(-k). 

(A8) 

When k = 0, we see W(k) is an antisymmetric matrix. 
As has been shown in Ref. |52j, the topological phase 
transition is determined by the Pfafhan of W(k) at k = 0. 

The Pfafhan of the antisymmetric matrix W(k) is di- 
rectly related to the topological invariance of the system. 
Dchnc 



k v\ — TJ ( ky 



= U 



(A3) 



where U is the SU(2) rotation matrix that takes the fol- 
lowing form, 



U 



cos(0) sin(0) 
-sin(0) cos(0) 



(A4) 



We can easily verify that UW = 1, therefore £kV = S,ka- 
On the other hand by defining a' z — a z , we hnd the new 
Pauli matrices a' x y z satisfy the standard commutation 
relations 



(A5) 



with s a bc the Levi-Civita symbol and 5 a b the Kronecker 
delta. The rotational symmetry plus the particle-hole 
symmetry ensure that the eigenvalues at each k should 
appear in pairs, thus in the topological superfluid regime, 
the lowest non-negative eigenvalue and the highest non- 
positive eigenvalue can become zero simultaneously and 
the Majorana Fermions are created in the vortex core of 
the topological superfluid. 

The above analysis cannot be applied to the case when 
h 7^ 0. Actually the in-plane magnetic field breaks the ro- 
tational symmetry of the effective Hamiltonian, because 
under the rotational symmetry in Eq. IA41 we have the 
following transformation, 



ha x — > h(cos(0)a x + sm(8)a y ). 



(A6) 



The breakdown of the rotational symmetry makes the su- 
perfluid asymmetric, rendering the unique FFLO vector 
Q, which is one of the major advantages of our proposal. 
The unique Q makes the detection of the FFLO vector 
more easily in realistic experiments, as discussed in the 
main text in sec. IIV1 



Q[ffrff(k)] = sign(Pf(W(k))). 



(A9) 



We see Q = ±1. Notice that Pf(A) 2 = Det(A) where A 
is any anti-symmetric matrix, and 

Pf(tf cff (k)A) 2 = Dct(tf cff (k)A) = Det(ff e ff(k)). (A10) 

Therefore Q changes sign only when Det(/f e g(k)) = 0. 
The deformation to the effective Hamiltonian generally 
does not change the value of Q, thus Q is topologically 
invariant. When k = and Q = (for BCS superfluid), 
we have 

Q[H eS (k = 0)] = sign^ 2 + T 2 - ( M 2 + A 2 )). (All) 

Thus for an in-plane Zeeman field, the system still has a 
topological phase transition, where the phase transition 
takes place at h 2 +T 2 = A 2 +fi 2 , at which point the band 
gap at k closes. We recover the well-known results when 
Zeeman field is set along the Z direction. Notice that 
in our model the order parameter and chemical potential 
should be determined self-consistently, where the Zee- 
man field generally destroys the superconductivity when 
h 2 + T 2 > A 2 , therefore the existence of the topological 
superfluid is not straightforward. 

We define the product of the eigenvalues, 



F = Y[E x =T>et(H efi ). 



(A12) 



For the system with the rotational symmetry (h = 0), we 
have 

F={h 2 + a 2 k 2 - A 2 - ^) 2 + 4a 2 fc 2 A 2 > 0. (A13) 

Obviously, F = happens at k = and T 2 = A 2 + ^ 2 . 
When h ^ 0, we find 

F= {h 2 -(A 2 +A-2\ r B))(h 2 -(A 2 +A+2^j). (A14) 
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Therefore F = gives 



h 2 = A 2 + A ± 2>/B. 



(A15) 



where A = 



a 2 (/c 2 - kl) - T 2 , and B = a 2 fc 2 (£ 2 - 
r 2 ) - a 2 fc 2 (A 2 + a 2 fc 2 ). When B < 0, F is always large 
than zero and the system is always gapped. We can easily 
verify that there is no real solution for A < and B > 0, 
therefore A > 0. Moreover, we find 

A 2 -4B = (r 2 -e 2 -a 2 (fc 2 +fc 2 )) 2 +4a 2 /fc 2 A 2 > 0. (A16) 

Hence A — 2\f~B > and h 2 > A 2 when the system has 
potential zero eigenvalues. Notice that for topological su- 
perfiuid, h 2 > A 2 + [i 2 is required. We find numerically 
that before reaching the regime h 2 > A 2 , the supcrfluid 
is always destroyed. Therefore for the BCS supcrfluid 
the system is always gapped and the topological phase 
transition can never happen. Notice that in our system 
the order parameter and chemical potential should be 
obtained self-consistently by solving the number and gap 
equations, therefore the existence of the topological su- 
pcrfluid is not straightforward. 

The sum of the eigenvalues has the following identity 



2k Q 



m 



(A17) 



Therefore for the BCS supcrfluid the sum of the eigenval- 
ues should be exactly zero, while for the FFLO supcrfluid 
the sum will never be equal zero, which means that for 
FFLO supcrfluid, the eigenvalues will never appear in 
pairs (E, —E). The above two results mark the impor- 
tant feature of the in-plane Zeeman field which differs 
dramatically from the Z direction Zeeman field. This 
is a direct consequence of rotational symmetry breaking 
due to the in-plane Zeeman field. 



Appendix B: Quantum Field Theory 

To facilitate the future investigation of the FFLO 
phase in ultracold atomic systems, we present the stan- 
dard field treatment of the FFLO phase in SO coupled 
degenerate Fermi gases. The same results in this ap- 
pendix have been obtained using the mean-field theory in 
the main text and these two methods give consistent de- 
scription of the FFLO physics. The quantum field theory 
has the potential of incorporating the fluctuation effect 
thus can be used in the future for studying the physics 
at finite temperature. 

The partition function at finite temperature reads as 



-S[ip,->pi 



(Bl) 



where 



Z = J V[1>,tf] 
Sty,^] = J drdv M^drM^)+H, (B2) 



where H = H a + VJ n t, with 

„2 

h = }2 ^(f-- /i+^ + ^Wi'. ( B3 ) 



CT,cr'=t4 



V so = a(p x a y - p y (T x ), and 



(B4) 



Here we assume a contact interaction, i.e., the two- 
particle interaction V(x, y) = gS(x — y) with g as the 
renormalized interaction strength. The FFLO order pa- 
rameter is 



<7<^(x)V t (x)) = Ae^-, 



(B5) 



where Q is the total momentum of the Cooper pairs, 
and A is a spatially independent constant. We make the 
following transformation, 



(B6) 



to remove the extra phase factor in the real space order 
parameter. Here we should emphasize that, in principle, 
we can assume ip a — > Vve^" x , where Q-f + Qj. = Q to 
make the pairing position independent. All these differ- 
ent choices of the Qi are equivalent class albeit the effec- 
tive Hamiltonian may be totally different. This property 
can be seen very clearly using the quantum field theory, 
but cannot be seen using the mean field developed in the 
main text. This flexibility may offer some advantage in 
dealing with some particular problems, however for the 
model considered in this work, we find Q^- = gives the 
simplest expression for the effective Hamiltonian, which 
is exactly the same as what we have obtained in the main 
text. For the single particle Hamiltonian, we find, 

<PtH aa , (p) W -+ (p)e iQ-x/2^ _ (B7) 

Thus we have the following transformation, 

ffaa'(p) -+ e- jQ - x/2 iW(pK Q - x/2 



# CTCT '(p + Q/2). 



(B8) 



The action without many-body interaction therefore 
reads as 

SobP,^} = f drdv ^d T Vv + V4^(P+§)?/V, 

<r=T4 

(B9) 

which means that we just need to replace p with p + Q/2, 
and, — p with — p + Q/2 so as to construct the effective 
Hamiltonian. Using the standard Hubbard-Stratanovich 
transformation, 

cxp(-g J dxdT^JvtytVa) = J V i A ] x 

></ 



exp( / drdx- 



AA* 



AVM + A*V^), (BIO) 
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the Fermi field now becomes quadratic thus can be inte- 
grated out, leaving only the bosonic field A and its con- 
jugate. The partition function Z — J T>Ae~ Sctt , where 
the effective action reads as, 



S e s = / drdr 
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— InDet/JG- 1 

2(3 



Tr(H), (Bll) 



where (3 = 1/T, Tr{H) = £ ff H ™> and G" 1 = d T + H eS> 
with the effective Hamiltonian 



-Heff = 



/ /C(p) A/ 2x2 
\A^I 2x2 -a y IC*(-p)a y 



Here 



JC(p) = 



i?tt(P ' 
H it(P 



Q/2) 
Q/2) 



#tl(p- 

H U (P 



Q/2) 
Q/2) 



(B12) 



(B13) 



The basis is arranged as (ipf, ip±, i/j], — ^>j). The above 
Hamiltonian is also identical to the effective Hamiltonian 
in Eq. [8]in the main text. For the particle-hole symmetry 
of the above Hamiltonian, see the discussion in the main 
text. 

For the saddle point approximation, we assume A = 
Ao+77, where Ao is spatially independent, while r\ — r/(x) 
corresponds to the fluctuation field. The contribution of 
the fluctuation field is not considered here, thus to the 
leading order, we have the effective action 



So 



|Ao| 



k,A. 



(B14) 



where ioo = (2n + l)irT is the Matsubara frequency. E\ 
is the eigenvalue of the effective Hamiltonian H c ff. Sum- 
ming over the Matsubara frequency using the standard 
method and the relationship between effective action and 
thermodynamical potential Vl = S c s/ (3 = S c gT, we have 



fi = — - 



9 2/3 



(B15) 



k.A 



At zero temperature, we have 



lim T\n(l + e- E / T ) = ° E " °' 
t^o+ v ; I -E E < 0. 



-EQ(—E), (B16) 



where Q(x) is the Heaviside step function. Thus we have 



CI = - 



|Ao| 2 , 1 
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YEx&(-E x ) + Y&- (B17) 



k,A 



We hereby recover the results from the mean field theory 
in Eq. [TTJ 



* Email: skylark.gong@gmail.com 
' Email: xbz@ustc.edu.cn 

* Email: chuanwei.zhang@utdallas.edu 

1 J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Phys. Rev. 
106, 162 (1957). 

2 P. Fulde, and R. A. Ferrell, Phys. Rev. 135, A550 (1964). 

3 A. I. Larkin, and Yu. N. Ovchinnikov, Zh. Eksp. Teor. Fiz. 
47, 1136 (1964). 

4 A. I. Larkin, and Yu. N. Ovchinnikov, Sov. Phys. JETP 
20, 762 (1965). 

5 A. Buzdin, Y. Matsuda, and T. Shibauchi, Europhys. 
Lett., 80, 67004 (2007). 

6 M. D. Croitoru, M. Houzet and A. I. Buzdin, Phys. Rev. 
Lett. 108, 207005 (2012). 

7 Y. Matsuda, and H. Shimahara, J. Phys. Soc. Jpn. 76, 
051005 (2007). 

8 K. Gloos, R. Modler, H. Schimanski, C. D. Bredl, C. 
Geibel, F. Steglich, A. I. Buzdin, N. Sato, and T. Ko- 
matsubara, Phys. Rev. Lett. 70, 501 (1993). 

9 A. Bianchi, R. Movshovich, C. Capan, P. G. Pagliuso, and 
J. L. Sarrao, Phys. Rev. Lett. 91, 187004 (2003). 

J. Singleton, J. A. Symington, M-S Nam, A. Ardavan, M. 
Kurmoo, and P. Day, J. Phys.: Condens. Matter 12, L641 
(2000). 



11 R. Lortz, Y. Wang, A. Demuer, P. H. M. Bottger, B. Bergk, 
G. Zwicknagl, Y. Nakazawa, and J. Wosnitza, Phys. Rev. 
Lett. 99, 187002 (2007). 

12 R. Casalbuoni and G. Narduli, Rev. Mod. Phys. 76, 263 
(2004). 

13 M. G. Alford, K. Rajagopal, T. Schaefer, A. Schmitt, Rev. 
Mod. Phys. 80, 1455 (2008). 

14 M. Houzet and A. Buzdin, Phys. Rev. B 63, 184521 (2001). 

15 L. G. Alsamzov, Zh. Eksp. Teor. Fiz. 55, 1477 (1968) [Sov. 
Phys. JETP 28, 773 (1969)]. 

16 A. V. Balatsky, I. Vekhter, J.-X. Zhu, Rev. Mod. Phys. 78, 
373 (2006). 

17 M. W. Zwierlein, C. A. Stan, C. H. Schunck, S. M. F. 
Raupach, A. J. Kerman, and W. Ketterle, Phys. Rev. Lett. 
92, 120403 (2004). 

18 C. Chin, R. Grimm, P. Julienne, and E. Tiesinga, Rev. 
Mod. Phys. 82, 1225 (2010). 

19 S. Giorgini, L. P. Pitaevskii, and S. Stringari, Rev. Mod. 
Phys. 80, 1215 (2008). 

20 P. Horak, J. Y. Courtois, G. Grynberg, Phys. Rev. A 58, 
3953 (1998). 

21 L. Guidoni, C. Triche, P. Verkerk, and G. Grynberg, Phys. 
Rev. Lett. 79, 3363 (1997). 

22 L. Sanchez-Palencia and M. Lewenstein, Nature Physics 6, 



13 



87 (2010). 

23 M.W. Zwierlein, J. R. Abo-Shaeer, A. Schirotzek, C.H. 
Schunck, and W. Ketterle, Nature 435, 1047 (2005). 

24 D. E. Sheehy and L. Radzihovsky, Phys. Rev. Lett. 96, 
060401 (2006). 

25 Y. L. Loh and N. Trivedi, Phys. Rev. Lett. 104, 165302 
(2010). 

26 H. Hu and X.-J. Liu, Phys. Rev. A 73, 051603(R) (2006). 

27 M. W. Zwierlein, A. Schirotzek, C. H. Schunck, and W 
Ketterle, Science 311, 492 (2006). 

28 M. W. Zwierlein, C. H. Schunck, A. Schirotzek, and W 
Ketterle, Nature 442, 54 (2006). 

29 G. B. Partridge, W. Li, R. I. Kamar, Y.-A. Liao and R. G 
Hulet, Science, 311, 503 (2006). 

30 W. V. Liu, F. Wilczek, Phys. Rev. Lett. 90, 047002 (2003) 

31 R. Winkler, Spin-Orbit Coupling Effects in Two- 
Dimensional Electron and Hole Systems, Springer Tracts in 
Modern Physics Vol. 191 (Springer- Verlag, Berlin, 2003). 

32 J. Ruseckas, G. Juzeliunal, P. ohberg, and M. Fleis- 
chhauer, Phys. Rev. Lett. 95, 010404 (2005). 

33 X.-J. Liu, M. F. Borunda, X. Liu, and J. Sinova, Phys. 
Rev. Lett. 102, 046402 (2009). 

34 X.-J. Liu, X. Liu, L. C. Kwek, and C. H. Oh, Phys. Rev. 
Lett. 98, 026602 (2007). 

35 C. Zhang, Phys. Rev. A 82, 021607 (2010). 

36 D. L. Campbell, G. Juzeliunas, and I. B. Spielman, Phys. 
Rev. A 84, 025602 (2011). 

37 J. D. Sau, R. Sensarma, S. Powell, I. B. Spielman, and S. 
Das Sarma, Phys. Rev. B 83, 140510(R) (2011). 

38 Y.-J. Lin, K. Jimenez-Garcia, and I. B. Spielman, Nature 
471, 83 (2011). 

39 Y.-J. Lin, R. L. Compton, A. R. Perry, W. D. Phillips, J. 
V. Porto, and I. B. Spielman, Phys. Rev. Lett. 102, 130401 
(2009). 

40 J.-Y. Zhang et ai, Phys. Rev. Lett. 109, 115301 (2012). 

41 C. Qu, C. Hammar, M. Gong, C. Zhang, P. Engels, in 
preparation. 

42 P. Wang, Z.-Q. Yu, Z. Fu, J. Miao, L. Huang, S. Chai, H. 
Zhai, J. Zhang, Phys. Rev. Lett. 109, 095301 (2012). 

43 L. W. Cheuk, A. T. Sommer, Z. Hadzibabic, T. Yefsah, 
W. S. Bakr, M. W. Zwierlein, Phys. Rev. Lett. 109, 095302 
(2012). 

44 L. He, M. Jin and P. Zhuang, Phys. Rev. B 74, 024516 
(2006). 



45 G. J. Conduit, P. H. Conlon, B. D. Simons, Phys. Rev. A 
77, 053617 (2008). 

46 Z. Zheng, M. Gong, X. Zou, C. Zhang and G.-C. Guo, 
arXiv: 1208:2029. 

47 J. Dalibard, F. Gerbier, G. Juzeliunas, and P. Ohberg, 
Rev. Mod. Phys. 83, 1523 (2011). 

48 V. Barzykin and L. Gor'kov, Phys. Rev. Lett. 89, 227002 
(2002). 

49 K. Michaeli, A. C. Pooter and P. A Lee, Phys. Rev. Lett. 
118, 117003 (2012). 

50 K. Martiyanov, V. Makhalov, and A. Turlapov, Phys. Rev. 
Lett. 105, 030404 (2010). 

51 A. P. Schyder, S. Ryu, A. Furusaki, and A. W. W. Ludwig, 
Phys. Rev. B 78, 195125 (2008). 

52 P. Ghosh, J. D. Sau, S. Tewari, and S. Das Sarma, Phys. 
Rev. 82, 184525 (2010). 

53 C. Zhang, S. Tewari, R. M. Lutchyn, and S. Das Sarma, 
Phys. Rev. Lett. 101, 160401 (2008). 

54 M. Gong, S. Tewari and C. Zhang, Phys. Rev. Lett. 107, 
195303 (2011). 

55 M. Gong, G. Chen, S. Jia, and C. Zhang, Phys. Rev. Lett. 
109, 105302 (2012). 

56 J. Zhou, W. Zhang and W. Yi, Phys. Rev. A 84, 063603 
(2011). 

57 M. Iskin and A. L. Subasi, Phys. Rev. Lett. 107, 050402 

(2011) . 

58 X.-J. Liu, L. Jiang, H. Pu, and H. Hu, Phys. Rev. A 85, 
021603(R) (2012). 

59 M. Iskin, Phys. Rev. A. 85, 013622 (2012). 

60 G. Chen, M. Gong, and C. Zhang, Phys. Rev. A 85, 013601 

(2012) . 

61 C. -H. Pao, S.-T. Wu and S. -K. Yip, Phys. Rev. B 73, 
132506 (2006). 

62 E. Altman, E. Dernier, and M. D. Lukin, Phys. Rev. A 70, 
013603 (2004). 

63 T. K. Koponen, T. Paananen, J-P Martikainen, and P. 
T6rma, Phys. Rev. Lett. 99, 120403 (2007). 

64 T. K. Koponen, T. Paananen, J-P Martikainen, M. r. 
Bakhtiari, and P. Torma, New. J. Phys. 10, 045014 (2008). 

65 J. Steinhauer, N. Katz, R. Ozeri, N. Davidson, C. Tozzo, 
F. Dalfovo, Phys. Rev. Lett. 90, 060404 (2003). 

66 L.-M. Duan, Phys. Rev. Lett. 96, 103201 (2006). 



